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UN-REDUCTION IN FIELD THEORY, WITH APPLICATIONS 


A. ARNAUDON, M. CASTRILLON LOPEZ, AND D.D. HOLM 


Abstract. The un-reduction procedure introduced previously in the context of Mechan¬ 
ics is extended to covariant Field Theory. The new covariant un-reduction procedure is 
applied to the problem of shape matching of images which depend on more than one 
independent variable (for instance, time and an additional labelling parameter). Other 
possibilities are also explored: non-linear u-models and the hyperbolic flows of curves. 


1. Introduction 

Symmetry (i.e., invariance under a Lie group of transformations) greatly facilitates the 
study of variational problems, both for the construction of explicit solutions of the vari¬ 
ational equations and for their qualitative analysis. A rich variety of information arises 
from Lie symmetry of variational problems, especially when they are formulated geomet¬ 
rically. For example, a vast, interesting literature exists on the topic of reduction by 
symmetry. In reduction by symmetry, the idea is to take advantage of the group of sym¬ 
metry transformations to reduce the dimension of the configuration and phase spaces of 
the variational problem, thereby making the problem easier to handle. When performing 
such a reduction, one must also provide a method of reconstructing the solutions of the 
original, unreduced, variational problem from solutions of the reduced problem, which 
sometimes requires additional compatibility conditions. 

Surprisingly, there are nice instances where this procedure can be used backwards. For 
example, suppose a variational problem looks complicated, but it may be recognised as 
the reduction by a certain group of symmetries of a variational problem formulated in 
a bigger space. Although the dimension of the corresponding un-reduced configuration 
space may be larger, the equations or the space itself may be simpler. Furthermore, the 
existence of the groups of symmetries may shed light on the nature of the initial equations. 
In this situation, one should notice that reduction by symmetries changes the structure 
of the equations. For example, in the Lagrange-Poincare reduction procedure (when the 
conhguration space is a manifold Q on which a Lie symmetry group G acts properly, see 
[5], [7], [14]), the reduced variational equations split into two different types. The first type 
is an Euler-Lagrange operator coupled with a gyroscopic term (the curvature of a chosen 
connection A in the bundle Q —)■ Q/G). The second type is a conservation law. In order 
to have a free variational problem in the reduced space, one needs to introduce forces into 
the un-reduced principle so that the equations will decouple. The choice of this force can 
be made by splitting the Lagrangian into horizontal and vertical parts with respect to the 
connection A. This is the un-reduction construction given in [4] for variational problems 
of a particle (Mechanics) and generalized in this article to a covariant field theoretical 
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setting. In particular, we also explore the topological situations which arise when the 
parameter manifold is not longer simply connected. 

The main motivation of [4] was shape matching: given two plane shapes 81,82 G 
Sh(M^), understood as closed curves in one seeks the optimal path of shapes joining 
81 and 82 - This problem is analysed in [9], [24] and references therein. The space Sh(M^) 
is a complicated inhnite dimensional manifold. However, we have Sh(M^) = Q/G, where 
G = Diff’'"(S'^), and Q is the space Emb^(S'^, M^) of positively embedded parametrizations 
of the circle in the plane, which is a much easier functional space than the unparame- 
terised planar curves in Sh(M^). By means of conveniently chosen forces, one may use 
un-reduction to lift the problem of shape matching to Emb’''(S'^,M^). In this article, 
this situation becomes richer. In particular, we can study matching of shapes depending 
on, say, two independent variables. A primary case is where the shapes depend on time 
(time evolution) and another parameter (space evolution) labelling a set of subjects in 
a research study. This so-called spatiotemporal analysis of shapes is a recent and active 
held of research. For details, the reader may consult [11], [15], [26]. In spatiotemporal 
shape analysis, there are two main approaches. These are the time-specihc and subject- 
specihc approaches, indicating the variable which parameterises the evolution in shape 
comparisons; either for a certain subject at a sequence of times, or for a sequence of 
subjects at a certain time. This spatio-temporal construction is illustrated in Fig. 1. 
Note that the x and t variables have interchangeable meanings. A more complex con¬ 
struction is found in [ 11 ] where the authors build a subject-specihc approach together 
with a time-reparametrization, with interesting applications to the compared evolution of 
Homo Sapiens Neanderihalensis and Homo Sapiens Sapiens, or bonobos and apes. The 
methodology is meant to couple with statistical analysis. The conhguration space of this 
approach is Diff(M^) together with the time reparametrization in Diff(M). 

The un-reduction procedure in Emb’'"(S'^,M^) that we propose here provides simulta¬ 
neous evolution of both types in a single system of partial differential equations. We 
expect this combined evolution to provide more accurate and versatile information for 
the problem of spatiotemporal curve matching. Furthermore, we introduce a certain con¬ 
venient Riemannian metric in the space of embeddings depending on derivatives of the 
curve (a Sobolev metric) which seems to be appropriate for the evolution in Sh(M^) and 
Emb’''(S'^,M^) (see [2], [3]), spaces which possess some other natural but pathological met¬ 
rics. In [1], the authors further investigated this approach with a simple numerical test 
in the classical mechanical setting, but more work is needed to obtain a reliable scheme. 

Because the theory is quite general, the range of potential applications is wide. Apart 
from the motivation of curve matching, we point out two other completely different areas 
of mathematical physics where covariant un-reduction is hidden. For example, cx-models in 
homogeneous spaces G/H may be written as an un-reduction problem in G. Interestingly, 
we may sometimes combine un-reduction with Euler-Poincare reduction to the Lie algebra 
0 to get a new set of equations. These equations are already in the literature, but we 
incorporate full geometric meaning to them with this concatenation of un-reduction and 
reduction, a situation intimately related with dual pairs (see [16]). Finally, covariant un¬ 
reduction is also applied to hyperbolic curve evolution, a baby geometric construction of 
other more sophisticated geometric flow equations. 


UN-REDUCTION IN FIELD THEORY, WITH APPLICATIONS 


3 



Figure 1. This diagram illustrates the spatio-temporal deformation of 
curves in Sh(S'^,]R^) that is considered in this work. The combination of 
spatial and temporal deformations, where the precise meaning of space and 
time has to be defined depending on the context, allows for a simultane¬ 
ous deformation of a curve along two parameters. The solution is then a 
function of {x, t) which minimises an given energy functional. In the sim¬ 
plest case of quadratic energy functional, the solution is known as being a 
harmonic map. 

Plan and Main Contents of the Paper. Section 2 reviews the basic concepts of 
covariant Lagrange-Poincare reduction, before formulating the main result of the paper, 
which is the Un-reduction Theorem 3.1, in Section 3. Section 4 provides examples of 
explicit applications of the Un-reduction Theorem for (i) curve matching in the plane; 
(ii) nonlinear sigma model and (hi) hyperbolic curve evolution. Each of these examples 
demonstrates the method of un-reduction and illustrates different ways to take advantage 
of the geometry of the reduced space. 

2. Covariant Lagrange-Poincare reduction 

The main result of the paper will be formulated as Theorem 3.1 in the next section. 
This section first reviews the basic concepts of covariant Lagrange-Poincare reduction. 
The version of this reduction in Mechanics takes place when a Lie group of symmetries 
G acts properly on the configuration manifold Q of the variational problem under study 
(for example see [7]). In the field theoretical setting, the group of symmetries acts on a 
fibre bundle tt : U —)■ iV by vertical diffeomorphisms, that is, actions such that 7^{y ■ g) = 
n{y),yy & E, g ^ G. We refer the reader to [5] and [14] for the exposition of the theoretical 
framework of this procedure. For our purposes, in this article we have adapted these 
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results as follows. On one hand we just consider trivial bundles Q x N ^ N, so that 
the dynamical objects of interest are mappings from N to Q and the problem is defined 
by a hrst order Lagrangian defined in the hrst jet space J^{N,Q). This simplification is 
mainly done for convenience in the applications, though the theoretical core of this work 
can be done in full generalities. On the other hand, we need to incorporate forces to our 
scheme, which will induce new terms in the equations in a straightforward manner. 

2.1. Background material. For the standard notions on bundles and connections, the 
reader can go, for example, to [19]; and to [18] for the basic definitions on geometric 
variational calculus in bundles and field theories. 

Let 7 r:Q—)-Q/G = Sbea G-principal bundle where the action Rg : Q ^ Q, g E G, is 
assumed to be on the right. Recall that a principal connection .4. is a 0 -valued 1-form in Q 
such that the equivariance property R^A = Ad^-i oA holds, and 41 (^q) = for any ^ G 0 , 
where is the infinitesimal generator of the action, i.e., ^Q^q) := d/de\^=oRexp{e^}{q) ■ This 
definition is equivalent to a choice of G-invariant splitting of the tangent bundle TQ into 
horizontal and vertical parts 

TgQ = HgQ © VgQ, 

for q E Q, where VqQ = ^ G 0 } and HgQ = ker4.. We denote by : TQ —HQ 

and p" ■. TQ ^ VQ the induced projections. The curvature of A is defined to be the 
0 -valued two form B = dA + [41,4.] and satisfies the equivariance property {Rg)*B = 
Adg-i o B. One can also define a 2 -form in S, but taking values in the adjoint bundle 
0 = (Q X q)/G as 

B{up, Wp) = [q,B{u^g, mJ)]g, Up, Wp E TpT, 

where Ug stands for the unique tangent vector (the horizontal lift of Ug with respect to A) 
in HgQ such that Tn^Ug) = Up. The definition does not depend on g G 7 r“^(p) because of 
the equivariant behaviour of the curvature. 

Let N be an oriented manifold endowed with a volume form v and consider a Lagrangian 
function L : J^{N,Q) -E M dehned in the 1-jet space of mappings s : N ^ Q. As 
the jet space R{N,Q) can be naturally identified with T*N © TQ, we will use both 
representations of this space in the following. We assume that L is invariant with respect 
to the lifted action of G in J^{N, Q), defined as 

:= jliRg o s) 

for g G G and any (local) mapping s. We can thus drop L to the quotient to obtain a 
reduced Lagrangian function 

i : J\N, Q)/G ~ T*N © {TQ)/G M. 

If we £x a principal connection A of the bundle Q ^ Q/G, we have a diffeomorphism 

(TQ)/G-^TS ©0 
(K]g) ^ {T7r{vg), [q,A{vq)]G), 
so that the reduced phase space decomposes as 

{J\N, Q))/G = T*N © {TQ)/G = T*N © (TS © 0 ) ^ J\N, S) © {T*N © 0 ), 
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SO that the reduced Lagrangian can then be written as 

£: J\N,E)(£){T*N^q) 

In the following sections, we will work with variational principles including a force 
term, that is, a map F : J^{N, Q) —)■ T*Q. The connection A splits the cotangent bundle 
T*Q = V*Q®H*Q and we can consider the decomposition F = F^+F"" where F^ = p^oF 
and F^ = p" o F with p^ and p^ denoting the projections of V*Q and H*Q respectively. 
We use the same notation as for the projection of the tangent bundle as no confusion can 
occur. If in addition F is G-equivariant with respect to the action of G in both the source 
and target spaces, we can drop F^ and to J^{N, Q)/G as 

f : J\N, S) © {T*N © 0 ) ^ T*E and /" : J\N, S) © {T*N © 0 ) ^ r • 

Note that for we use H*Q/G ~ T*E, and for we have the isomorphism VQ/G ~ 0 
given by ^ [g,'^]G- 

Finally, we recall the dehnition of the canonical momentum map for the natural lift 
action of G on T*Q 

J;T*Q ^ 0 * 

(J(«g),OgX0* = {piqAQ)TQxTQ* 

where Ug G T*Q, ^ G 0 , and G TQ. We can extend J to a map 

i -TN ®T*Q ^TN (2.1) 

trivially in the factor TN. We note that, if we identify TN A'^~^T*N, n = dimA^, by 
means of a fixed volume form v, the map J : TA^ © T*Q -A TN ® q* is the covariant 
momentum map in field theories (cf. [18] [Proposition 4.4]). 

2.2. Lagrange-Poincare reduction. In the sequel, we assume that N is compact. If 
N is not compact, the domain of variations of the maps s : N ^ Q will be assumed to be 
compactly supported. We project the variational principle defined for L from J^(A^, Q) to 
its quotient J^(A^, Q)/G. for £ : J^(A^, S) x (T*N © 0 ) —)■ M. Critical solutions are maps 
a : N ^ T*N © 0 which, moreover, project to maps p\ N^T = Q/Gasp = 'K^oa 
according to the diagram 

T*N 

/ ( 2 . 2 ) 

AT -4 S 

where TTg : T*N ©0 —)■ S is the projection of the adjoint bundle forgetting the T*N factor. 
The free variations of the initial problem provide a family of constrained variations that 
define a new type of variational equations. They are called Lagrange-Poincare equations 
(see [5], [14]). The next theorem gives the Lagrange-Poincare reduction with forces F 
which is obviously the one in the literature when F = 0 . 

Theorem 2.1 (Covariant Lagrange-Poincare reduction with forces). Let tt : Q —)■ 
Q/G = T, be a principal G-bundle, A be a principal connection on it and N be a compact 
manifold oriented by a volume form v. Given a map s : N ^ Q, let a \ N ^ T*N ® q be 
defined as 

a{x){u;) = [s{x),A{Ts ■ (a;))]G, 
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with u G TxN,x G N; and let p : N ^ T., p{x) = [s(a;)]G = TTg o a. We consider a G- 
invariant Lagrangian L : J^{N,Q) —)■ M and a G-eguivariant force F : J^{N^Q) —)■ T*Q. 
Then the following points are eguivalent: 


(1) s is a critical mapping of the variational principle 


6 [ L{s,jh)\+ [ {F{s,jh),6s)\ = 0 (2.3) 

Jn Jn 

with free variations 6s. 

(2) The Euler-Lagrange form of L satisfies the relation 

^£(Lv)(/s)=F. 

(3) a : iV —)■ T*N^ q is a critical mapping of the variational principle 


Jn 


[ {f{j^PW),Sp)^ + 

Jn 


{rifpw),v)^ = 0 , 


>N 


for variations of the form 6a = V-^p — [a^p] + B{6p,Tp) G 0 , where 6p G TpE is 
a free variation of p and p is a free section of g ^ T,. 

(4) a satisfies the Lagrange-Poincare eguations 






6i 


6i 




6a' 
6i 




6 a 


6 a 


(2.4) 


where SCp{£N) : J^(iV, S) —T*S is the Euler-Lagrange form of £ with respect to 
the variable p only and div^ stands for the covariant divergence operator defined 
by the connection A. 


Remark 2.2. Given a solution of the Lagrange-Poincare eguations (2.4), the reconstruc¬ 
tion of a solution of the initial variational problem reguires a compatibility condition. 
Given the map a : N ^ T*N 0 g and the induced map p : N ^ T,, we consider the pull¬ 
back principal bundle p*Q —)■ N and the pull-back of the connection p*A. Recall that the 
space of connections is an affine space modelled over the vector space ofg-valued 1-forms 
in the base manifold. We can thus consider the new connection A"^ = p*A-l-a. Then, the 
compatibility condition is 

Curv(^^) = 0. (2.5) 

Indeed, if this condition is satisfied, and the manifold N is simply connected (see ^3.4 
below for some topological issues), then the solutions s : N ^ Q are the integral leaves or 
sections of that connection. See [5, 14] for more details. 


3. The covariant un-reduction scheme 

We are now almost ready to describe the un-reduction scheme for Field Theories. As 
in the case of Mechanics (see [4]), this construction requires that the Lagrangian is de¬ 
composed into horizontal and vertical parts with respect to the chosen connection A. 
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3.1. Vertical and horizontal Lagrangians. We first give an expanded expression of 
the Euler-Lagrange form SC{L) : J‘^{N,Q) —)■ T*Q for an arbitrary Lagrangian L : 
J^{N, Q) —)■ M once a linear connection V in Q has been hxed. For that, we consider the 
horizontal lift v ^ v from TQ to T{T*N ® TQ) with respect to V (the lift is done in the 
TQ part only and is trivial in the T*N factor). Then we dehne ^ : J^{N, Q) —)■ T*Q as 


:=dT(j»-(5s, 

\ / TQxT*Q 

for any 5s G TqQ, q = s{x). On the other hand, we dehne the vertical derivative : 
J\N,Q) TN®T*Q as 




L{jls + eu), 


e=0 


for any u E T*N <S) TgQ, q = s{x). The Euler-Lagrange form is thus 


SC{L){fs) = (3.1) 

where div^’'" stands for the divergence operator dehned by the volume form v and the 
affine connection V. It acts on valued vector helds in N (note that along the map 
j^s, dL/dj^s is precisely a section of TN ® s*T*Q —)■ N) and it is dehned as the only 
operator such that 

div'^ (T’,X) = {|div^’^A’,X^ + (V, VX) 


for any vector held X eTN ® T*Q and any section vector held X in TQ. 

We now assume that the Lagrangian L : J^(iV, Q) = T*N®TQ —)■ M can be decomposed 
as L = L^ + L^ with 


L^{u v) = L^{u p’^{v)) and {u v) = {u p"(v)) 

for any u ^ v E T*N ® TQ, with respect to the connection A. Furthermore, as TQ = 
HQ (BVQ, we have 

: T*N ^HQ^R and H : T*N ®VQ^R. 


Obviously, the G invariance of L and A extends to the G-invariance of U and so that 
they drop to the quotient as 

t : J^{N,T)=T*N ®TT^R and G:T*N®~q^R, 
to form the reduce Lagrangian i{j^p, <j) = + G{p, a). It is easy to see that 

5£ _ 5£^ 5£ _ 5G 

5j^p 5j^p 5a 5a 


We then consider that the linear connection V in Q is invariant under the action of G so 
that it projects to a linear connection V in S = Q/G by the condition VxY = T^*{yx^Y^)- 
In addition, the connection A induces a connection in the associated bundle 0 —)■ S. With 
respect to these connections we can compute 

W£ V£^ vr 
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and the Lagrange-Poincare equations (2.4) thus read 


div^’" 






5j^pJ 6p 

6i 


= /" + 


div^^ + ad;^ = / 




vr 

6p 


6a 


6a 



(3.2) 


The Lagrangian splitting is crucial in this methods and allows the appearance of the 
standard Euler-Lagrange equations for & in the left hand side of the hrst equation. The 
second important ingredient is the force term which will allow us to exactly obtain the 
Euler-Lagrange equations by cancelling the right hand side of the same equation. 


3.2. The un-reduction theorem. We are now ready to state the central theorem of 
the un-reduction method using the held theoretical context described above. 


Theorem 3.1. Let N be a smooth manifold oriented by a volume form v and tt : Q —)■ S 
be a G-principal bundle equipped with a principal connection A. Let I : = 

T*N (8 TS —)■ M 6e a first order Lagrangian. We consider a G-invariant Lagrangian 
L : Q) = T*N ® TQ —)■ M such that L = L^ + where L^ o p^ = L^ is uniquely 

determined by I, L^ o p^ = L^ is arbitrary, and p^,p" are the projectors of the splitting 
TQ = HQ (BVQ induced by A. We also consider a G-equivariant force F : Q) —)■ 

T*Q such that F"" = p" o F is arbitrary and F^ = p^ o F is given by the condition 


vr , /6G . - 


(3.3) 


for its projection S) x {T*N 0 g) —)■ T*S. Then, the variational equations of 

the problem defined by L and F read 

£C{L^){jh) = 0 

where .4,* : g* —)■ V*Q is the dual of the connection form. Finally, critical solutions 
s : N ^ Q of (3.4) project to critical solutions p = [s]g of the Euler-Lagrange equations 
SC{l){fp) = 0. 



Proof. We follow the notations of the preceding sections. The variational principle of L 
and F is 


0 = 6 I L’^v+ 6 I LW+ (F\(5s)v+ / {F^,6s)v 

J N 


' N 


' N 


>N 




>N 


6P 


+ / (/*Xp>v+ / {F\A{Ss)) 


ve 

N \ 

V 


'V 


IN 


IN 


IN 


IN 


6L 


= 6 L{js)n+ ( —,i3(Tp,V) )v 


6a 


-,6p)n 


6t 


IN 


6a 


In 


{F\Ai6s))N. 
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From the expression of 6 a in Theorem 2.1 with ri{x) = (s(a;), ^(5 s))g we have that 




6 a + B{Tp, 6p) 


V = 




,V' 


A 



V. 


For any f : N ^ q, we recall that the covariant derivative is V'^(s {x),f{x))G = 
{s{x),df{x) + [A{j^s),f])G = {s{x),df{x))G + [a, (s(a;),/( x))g]. Now, for / = ^(5s), 
we have 




6t - \ f 16t \ 

— ,5cr + i3(Tp, V)y V = J ^ —,(s,d^(5s))Gy V 

,dA(S.)y=-j^ ^div'' J V, 


Finally note that, as L^{j^s) = the variation of the action defined by with 

respect to vertical variations of antomatically vanishes. The variational principle 
natnrally splits into vertical and horizontal part as equations (3.4). 


Solutions of the variational problem defined by M project to solutions of the problem 
dehned hj I = by Theorem 2.1. □ 

Remark 3.2. If we have iV = M, v = df, (that is, the case of classical Mechanics) we 
have div'" = d/dt and we recover the results and equations of [4]. 


The expression of the horizontal force defined by condition (3.3) is 


ph 


ds 



The variational principle on the un-reduced space of equation (2.3) is then defined using 
this particular force such that the reduced Lagrange-Poincare equations decouples. 


The hrst equation in (3.4) is the usual Euler-Lagrange equation for the horizontal 
Lagrangian. With respect to the second, we first note that we cannot exchange the 
position of A and div'" as the authors do in [4]. In fact, the divergence of A*3{6L'^/6j^s) 
would require an additional (linear) connection in Q. Moreover, as we mentioned in the 
definition (2.1) of J above, we have that J{6L'^/6j^s) is a covariant momentum map, so 
that div^3{6M/6j^s) is the expression of a conservation law with respect to the group of 
symmetries. The second equation in (3.4) equals this to the vertical part of the force. If 
one set F'^ = 0, the conservation law is complete, but sometimes it is interesting to keep 
this vertical force as it might be used to externally control the dynamic along the vertical 
space. 


3.3. Reconstruction and the surjectivity of the un-reduction scheme. The the¬ 
orem 3.1 above says that solutions of the un-reduced problem project to solutions of the 
Euler-Lagrange equations dehned by the Lagrangian 1. One may ask if this projection is 
exhaustive, that is, if every solution of the variational equations of Z is a projection of a 
solution of L. This question involves some topological constraints concerning N (see §3.4), 
but we hrst give an answer assuming that N is simply connected (or we just consider the 
question from a local point of view). From the Lagrange-Poincare reduction theorem, the 
variational equations dehned by L are equivalent to 

Xfv xpv 

= 0. div"*—+ ad;—= /•, 






10 


A. ARNAUDON, M. CASTRILLON LOPEZ, AND D.D. HOLM 


that is, they contain the Euler-Lagrange equations oi I = ^^ together with an additional 
set of equations which might restrict the solution of the hrst set. They are written in 
terms of the map a \ N ^ T*N ® q and p : N ^ T,. Recall that a determines p as 
p = TTg o cr (see diagram (2.2)). The key point is that the hrst reduced equation only 
involves p and its hrst jet j^p. Once we have a solution p and we may consider both 
the second reduced equation and the compatibility condition. They are now equations for 
maps a seen as sections of the bundle T*N ® p*Q —>■ N ^ which means that we “restrict” 
the vertical part of our construction to the hbers which sits only on the solution p on the 
base manifold. With the solution of these last equations, we can perform reconstruction 
to get a map s : N ^ Q such that p = [s]g- Roughly speaking, the reduced equations 
are uncoupled, so that p and a can be treated separately and the surjectivity of the un¬ 
reduction technique is guaranteed. The reason of this is the force term which exactly 
decouples these equations although it is not explicit in the un-reduced equations. 

3.4. Topological constraints and un-reduction. The topology of the manifold N 
may create interesting situations in the reconstruction and un-reduction frameworks. If 
N is not simply connected, the flatness of a connection, that is the compatibility condition 
(2.5), does not ensure the existence of global integral sections and the surjectivity of the 
projection s i—)■ p of solutions involves some other global considerations. 

An example of this situation is the following. Consider Q = and G = so that 
Q —)■ S = 5^ is the Hopf hbration. Choose the mechanical connection A in this bundle, 
that is, the connection such that HgS^ T VgS^ with respect to the standard Riemannian 
metric in S^. For the sake of simplicity we consider N = S^, that is, a problem of 
Mechanics with cyclic solutions where, in addition, the compatibility condition (2.5) is 
satisfied automatically. We denote 6 the coordinate of and we consider the G-invariant 
Lagrangian L : J^{N, S^) — )■ M, 

where s = ds/dO G Ts(e)S^, as well as its decomposition L = + U induced by the 

orthogonal splitting s{9) = s^{9) -f s'"{9) defined by the A. The adjoint bundle 0 —)■ 
is a trivial line bundle and the reduced phase space J^(iV, S) x {T*N ® g) becomes 
TS'^ X T*S^. We can write the reduced Lagrangian as ^ = G + G with 



where p : —)• S = S'^, p = dp/dO, and a = <,d6 with a map g = M. The reduced 

equations are 

vp = 0, = r- 

Solutions of the hrst equation are closed geodesics p in S^. Given one of these, the curves 
s{6) of the un-reduced problem will be in the restriction p*S^ of the Hopf hbration along 
p. This restriction is a torus and according to the reconstruction process seen in §2.2, 
the curve s{6) must be horizontal with respect to the connection A + <,d9. Under these 
circumstances, the curve s{6) need not be closed and in fact, the phase p G such that 
s(27r) — s(0) = p is precisely the holonomy of the connection A + ^dO along the curve p. 
The holonomy of A alone is vr (indeed, the connection A is not hat and the holonomy is 
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related with the Chern number of the Hopf bundle, see [19, Chapter XII]). Hence, besides 
conditions ^ = f" and ?(27r) = <^(0), for the closeness of c{6) we need <^{6) to satisfy 


- 27 r 


/ q{9)d9 = —TT, 

Jo 

so that we cancel the holonomy of A. Only very specihc functions may accomplish these 
conditions. For example, P{9) = cos(6*) gives <^(6*) = sin(6') — 1/2 as a possible solution. 
Other functions does not provide closed curves c{9). Furthermore, it is important to 
note that the constant value of the holonomy of the hxed connection A along geodesics p 
is unusual and other choices of hxed connections A will dehne a holonomy depending on 
p. In that case, the choice of f" will depend on the global curve p and will not be a local 
object. 

In other words, there are circumstances where one cannot recover all solution of the 
reduced problem from those of the un-reduced problem. It seems that the freedom in 
the choice of L'’ and, especially, might solve this issue but their specihc expression 
will depend on the solution p itself. We refer the reader to [23] and [25] to some related 
approaches to the problem or, for example, [28] for a similar situation to the example 
above in the context of isoholonomic problems and quantum computation. The situation 
for manifolds N of dimension greater than 1 is, of course, much more complicated. 


4. Applications 

4.1. Planar curve matching. We begin the application section with curve matching, 
he main motivation of this work, initiated by [4] and extended here to held theories. 

4.1.1. Geometric setting. Let Q = Emb’''(S'^,M^) be the manifold of positive oriented 

embeddings from to Elements in Q are maps c{9) G for 9 & and elements 
in the tangent space TcQ are pairs (c,n) with c G Emb^(S'^,M^) and u G a 

parametrized vector held along the curve c. Then 

TQ = Qx C°°{S\R^) 

and we can take a trivial linear connection V in Q. We consider an open domain iV C MxM 
with the Euclidean metric, coordinates {t, x) and volume form \ = dt A dx. Elements of 
the jet bundle J^{N, Q) ~ T*N (g) TQ are written as 

= ct{9){t,x)dt + Cx{9){t,x)dx, (4.1) 

that is, Ct and are the derivatives of a map c : N ^ Q along t and x respectively. 

We now consider the group G = Dih’''(S'^) of orientation preserving diheomorphisms 
of and its Lie algebra g = A(S'^) which consists of vector helds on S^. The group G 
acts on the right in Emb^(S'^,M^) as reparametrization of curves c and the reduced space 
is the space of shapes in 


^ Q _ Emb+(5^M2) 

''~G~ Dih+(^i) ' 


(4.2) 


The principal bundle Q —)■ S is endowed with a canonical principal connection A as 
follows. Given u G TcQ, we consider its tangent and normal decomposition 

u{9) = v{9)t{9) + h{9)n{9), 
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where (t, n) is the orthonormal Frenet frame along c and v{6), hid) scalar fnnctions along 
the cnrve. We clearly have that v{9)ti9) is a vector tangent to the orbits of G = Diff’''(S'^) 
so that vi9)ti9) G VcQ. We can thus define the horizontal part of u as the part h(6')n(6') 
and we have a decomposition TQ = HQ Q)VQ. 

The dehnition of a convenient Riemannian metric in Q = Emb’''(S'^,M^) invariant 
with respect to the action of G = Diff’''(S'^) is an interesting topic which has attired the 
attention of many research works (see, for example, [2] and [3] and the references therein). 
The natural metric 

9{u,v)=[ {u{9),v{9))dl, (4.3) 

Js^ 

with u,v E TcQ, and dl = \ce\d9 being the arc-length, is not very useful as it defines a 
zero geodesic distance in both Q and Q/G. The problem can be overcome in the shape 
space Q/G hy the metrics 

g{u,v)=[ (1 +AK,{9)‘^){u{9),v{9))dl, (4.4) 

Js^ 

with R > 0 and k the Frenet curvature of c. But this metric defines again a zero geodesic 
distance in Q along the hbers of the hbration Q ^ Q/G. A metric with a well dehned 
Riemannian distance in both Q and Q/G is obtained by adding higher order derivatives 
of u and n in a Sobolev-type expression as 

g{u,v) = [ {{u{9),v{9)) + A‘^{D0u{9),Dev{9))) dl, (4.5) 

Js^ 

where Dq = i^\de is the arc-length derivative. We can collect these three cases (as well 
as many others, see [2]) as 

9 v{u,v)=[ {u{9),Vv{9)), (4.6) 

Js^ 

for a convenient choice of a G-invariant self-adjoint pseudo-differential operator V which 
can depend on the curve and its derivatives. In particular, the operator for (4.4) if 
"P = 1 -|- Ak^, and for (4.5) we have V = 1 — A^Dg. One additional advantage of the 
operator associated to (4.5) is that it does not depend on the curve, whereas the operator 
for (4.4) depends on the curvature of the curve where it is evaluated. This represents a 
great simplihcation in the expression of the un-reduced equations. 

Remark 4.1. Even if the mechanical connection A in this context of space of embeddings 
is easy to visualise and compute, its structure is far from being trivial. The calculation of 
its holonomy and curvature is a whole subject on its own, which should be addressed in 
forthcoming works. 

4.1.2. Reduction and un-reduction. Elements of the shape space of plane curves S = 
Emb’''(S'^,]R^)/Diff“''(S'^) are denoted by p and the elements of the jet space AiN, S) = 
T*N ® TS are expressed as 

jlt,x)P = Pt(ty ^)dt + Px{t, x)dx. 

Furthermore, elements of T*N ® g are 

ait, x) = atit, x)dt + a^it, x)dx 
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where at{t,x),axit,x) belong to the adjoint bundle g —)■ S and can be understood as 
vector helds along a shape p G S and tangent to it. We consider the Diff'^(S'^)-invariant 
Lagrangian L : J^{N, Q) ~ T*N (g) TQ —)■ M 

= + (co;, Vc^)) dl (4.7) 

which can be decomposed a.s L = + L"" with respect to the connection A as 


L^Uix,tf) = \ [ 'Pht) + {h^, VK)) dl, 

L^’Hlx^tf) = ^ / Pvt) + (n^, Vv:,)) dl, 

where 

Cf = Vtt + and = vA + h^n. 


The un-reduction equations (3.4) are then computed in the proposition 4.2 below in 
the case when V is independent of the curve. 

Proposition 4.2. The un-reduced equations (3.4) for the hi-dimensional problem of pla¬ 
nar simple curves defined by the Lagrangian (4.7) and the metric (4.6) withV independent 
of the curve are 

dxPhx + dtPht = Do{hxPvx + htPvt) - kH 

dxPvx + dtPvt = ^ 

with the decomposition 

Cx = Vxt + hxiv, ct = vA + htu 

and for any choice of vertical force , where 

H =^{hxPhx + htVht). (4.9) 


Proof. The Euler-Lagrange equation contains two terms, the first is readily 

5L^ 

div-— = dtifPht) + dxiPhx). 

bAc 

Before computing the second term of the EL equation, we rewrite only the temporal part 
of the Lagrangian in order to simplify the calculation as 


L’^{c,Ac)\t = \ ■ n)n,P(ci ■ n)n)d/ = ^ j 


ce 


ce 


Ce 


CfJi^]P{cfJi^]\ce\de. 


Ce 


This Lagrangian being horizontal, we just need to consider variations of c that are hori¬ 
zontal with respect to A, that is, variations of the form 6 c = n,^, f G C°°{S^). With the 
identities Den = —nt, Jn = — t and Jt = n we compute 

f)fh\ r 

■ (nO = {ct ■ J (nOe) P (q ■ n) 


dO 


= fe{cf Jn)V {ct ■n)d9 + 
Js^ 


'SI 


f (ct ■ JDen) V (q ■ n) dl 


= / fDe [{ct ■ t) V {ct ■ n)] dl - fn (ct ■ n) P (q ■ n) dl. 


'si 


'si 
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Therefore the derivative of the Lagrangian is 

= DeihtVvt) - KhtVht. 
oc 

From the symmetry t ^ x, the Euler-Lagrange equation follows. □ 

Remark 4.3. The term kH can he interpreted as a penalty term for the deformation of 
most curved regions of the curve. The sign of this term would depend on the concavity or 
convexity of the curve at this point, and thus this force would try to prevent the curve to 
be deformed too fast in these regions. The eguation (4.8) also shows that the dynamics in 
{x, t) is governed by the coupling between ht and Vt reguired for the shape deformation to 
be independent of the reparametrisation. 

Remark 4.4. The un-reduced eguations with curvature weighted metric (4.4) can be com¬ 
puted directly from the variational principle, as in [4] . The eguation will have the same 
symmetry x -H- f but with more complicated terms. Because this metric is not very useful, 
the covariant un-reduced eguations will not be displayed here. 

4.2. Horizontal Lagrangians and cx-models. The freedom in the choice of forces and 
Lagrangians in Theorem 3.1 permits the trivial choice of = 0 and F"" = 0. From (3.3), 
the horizontal part of the force automatically vanishes. This simple situation appears 
when the un-reduced Lagrangian L is just the pull-back of the Lagrangian —)■ 

M with respect to the projection J^{N,Q) —)■ J^(A^, S), j^s h-)■ = i^P- A solution 

of the problem defined by L is any map s : N ^ Q such that p = [s]g is a solution for 
1. This means that there is a gauge degeneracy in the sense that, given a solution s and 
any map g : N ^ G, the map s = s ■ g is also a solution. 

Even though these trivial choices for F and U’ are not always convenient, there are 
some instances where they appear naturally. This is the case of cr-models in homogeneous 
spaces (see for example [12], [13], [17], [10]). Let Q = Ghe a. Lie group and FT be a closed 
subgroup such that the Lie algebra decomposes as 0 = m © f) for certain vector space m 
such that [P),m] Cm (that is, we have a reductive decomposition). We can right translate 
the decomposition m©f) = 0 = Tf,G to every TgG, g E G, thus obtaining a connection A for 
the principal bundle G ^ M over the homogeneous space T, = M = G/H. We consider 
the harmonic, or cr-model, problem on maps p : N ^ M defined by the Lagrangian 

1:J\N,M) M 

where the norm is taken with respect to a pseudo-Riemannian metric in N and a Rie- 
mannian metric in M. The lift L of / to J^{N,G) is 

L-.J\N,G) M 

3^9 ^ \\\p^{d9)f ^ 

where p^ : TG —?• HG is the horizontal projection defined by A and the norm is taken 
with respect to the metric of N and the lift of the metric in M to horizontal vectors in 
M. Theorem 3.1 can apply and solutions of the force-free problem defined by L project 
to the desired harmonic maps in M. 
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In the majority of the homogeneous spaces where relevant cx-models are dehned, the 
group G is endowed with a bi-invariant metric. In this case, the reductive decomposition is 
assumed to be m = ()■*■ and we have a metric in M by imposing the projection n : G ^ M 
to be an isometric submersion, that is, the metric in T^M is the same as the metric in 
HgG for any g with '^{g) = x. The group G left-acts on the coset space M by isometries. 
Hence, the Lagrangians I and L are both G invariant. This group of symmetries is too 
big for M to do reduction (in fact the orbit space is a single point), but we can perform 
covariant Euler-Poincare reduction for L. We then get a new reduced Lagrangian 

1: J\N,G)/G = T*N M 

where is the splitting dehned by the reductive decomposition. It is easy to see 

that the Euler-Poincare equations are 

div''<^m -F [<^0, <tm] = 0 

which, together with the suitable compatibility condition, can be used to get solutions of 
L that, afterwards, can be projected to S. This approach is found, for example, in [13], 
[10], [21]. The advantage of this un-reduction and reduction procedure relies on the fact 
that 0 is a simpler space (is a vector space) than either G and M. 

The situation can be even put in a more general framework as follows. Let L be a 
hrst order Lagrangian on a Lie group G as conhguration space, which is right invariant 
under the action of a subgroup H and left invariant under the group G itself. Suppose 
that we are interested in the induced variational problem in the homogeneous space G/H. 
The un-reduction and reduction procedure will give hrst a variational problem in G to 
hnally induce a problem in the Lie algebra 0 which, in general, is simpler. See [27] for a 
description of a similar situation in Mechanics (that is, = M). 

4.3. Hyperbolic curvature flow. The hyperbolic curvature flow of plane curves (see 
for example [22] or [29]) is the variational equation defined by the Lagrangian 

L : TEmh{S\R‘^) M, 



Note that this is not a geodesic variational principle of the metric (which provides 
null geodesic distances in both the curve and shape spaces) but a Lagrangian involving 
a kinetic and a potential term. Moreover, the Lagrangian L can be easily split into 
horizontal and vertical with respect to the connection A{ct)= {ct ■ t)t as 



where 


Ct = hn + vt. 

The Lagrangian L (and U’) is Diff(S'^)-invariant as its definition is geometric and 
does not depend on the parametrization of c but only on its shape. 
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One of the main featnres and applications of the hyperbolic flow (as well as of other 
geometric flows of curves) is the study of the evolution of the shapes of the curves under 
it. We can now suppose that we just want to study this evolution in the shape space 
Emb(S'^,M^)/Diff’'“(S'^). The natural Lagrangian in this situation becomes I = the 
projection of to this quotient space. In this context, the un-reduction technique applies 
and we have the last result of this paper. 

Proposition 4.5. The un-reduced equations for the system described above read 

dth = Deiyh) — — 1), 

dtv = 

Ct = hn + vt. 

In particular, if we choose F'" = Q and the initial tangent velocity to vanish (n(0) = Oj, 
then v{t) = 0 for all times and the velocity of h is proportional to the curvature n. 

Proof. Following the derivation of the un-reduced equation for curve matching, but in the 
classical case, one can prove this proposition as well. □ 

Remark 4.6. The equations of the Proposition above for F'" = 0 are the hyperbolic mean 
flow equations (see for example [22] for a good account of this flow). The usual approach 
in the literature works in Emb(S'^,M^) and then restrict oneself to the normal part of the 
flow. The approach here works with shapes in Emb(S'^M^)/Diff(S'^) so that the trivial 
choice of F'" = 0 gives directly the geometric equations. 
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